home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CPBFA.z / CPBFA
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCPPPPBBBBFFFFAAAA((((3333FFFF))))                                                            CCCCPPPPBBBBFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CPBFA   - CPBFA factors a complex Hermitian positive definite matrix
  10.      stored in band form.
  11.  
  12.      CPBFA is usually called by CPBCO, but it can be called directly with a
  13.      saving in time if  RCOND  is not needed.
  14.  
  15.  
  16. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  17.       SUBROUTINE CPBFA(ABD,LDA,N,M,INFO)
  18.  
  19. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  20.      On Entry
  21.  
  22.      AAAABBBBDDDD COMPLEX(LDA, N)
  23.         the matrix to be factored.  The columns of the upper
  24.         triangle are stored in the columns of ABD and the
  25.         diagonals of the upper triangle are stored in the
  26.         rows of ABD .  See the comments below for details.
  27.  
  28.      LLLLDDDDAAAA INTEGER
  29.         the leading dimension of the array  ABD .
  30.         LDA must be .GE. M + 1 .
  31.  
  32.      NNNN INTEGER
  33.         the order of the matrix  A .
  34.  
  35.      MMMM INTEGER
  36.         the number of diagonals above the main diagonal.
  37.         0 .LE. M .LT. N .  On Return
  38.  
  39.      AAAABBBBDDDD an upper triangular matrix  R , stored in band
  40.         form, so that  A = CTRANS(R)*R .
  41.  
  42.      IIIINNNNFFFFOOOO INTEGER
  43.         = 0  for normal return.
  44.         = K  if the leading minor of order  K  is not
  45.         positive definite.  Band Storage
  46.         If  A  is a Hermitian positive definite band matrix,
  47.         the following program segment will set up the input.
  48.         M = (band width above diagonal)
  49.         DO 20 J = 1, N
  50.         I1 = MAX0(1, J-M)
  51.         DO 10 I = I1, J
  52.         K = I-J+M+1
  53.         ABD(K,J) = A(I,J)
  54.         10    CONTINUE
  55.         20 CONTINUE LINPACK.  This version dated 08/14/78 .  Cleve Moler,
  56.      University of New Mexico, Argonne National Lab.  Subroutines and
  57.      Functions BLAS CDOTC Fortran AIMAG,CMPLX,CONJG,MAX0,REAL,SQRT
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.